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Abstract 



Thermodynamic fluctuation theory originated with Einstein who inverted 

the relation S* = fc^ In $7 to express the number of states in terms of entropy: 

I Q = exp{S/kB)- The theory's Gaussian approximation is discussed in most 

'^ statistical mechanics texts. I review work showing how to go beyond the 

Q Gaussian approximation by adding covariance, conservation, and consistency. 

O This generalization leads to a fundamentally new object: the thermodynamic 

Riemannian curvature scalar R, a thermodynamic invariant. I argue that \R\ 
cn is related to the correlation length and suggest that the sign of R corresponds 

J!C to whether the interparticle interactions are effectively attractive or repulsive. 

^. 1 INTRODUCTION 

Q Thermodynamic properties follow from particle properties via statistical mechanics. 

^-H But thermodynamics is more than an application of statistical mechanics. Thermo- 

>■ dynamics can sometimes yield microscopic information by means of the fluctuation 

k> theory as originated by Einstein in 1904 pLj. 

5-1 What does thermodynamic fluctuation theory tell us? In its well-known Gaussian 

approximation it yields the probability of finding a system in some thermodynamic 
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state. Most texts do not go further than this Gaussian approximation. However, 
here I argue that going further using ideas from Riemannian geometry adds a fun- 
damentally new element to the picture: the thermodynamic Riemannian curvature 
R, which yields information about interparticle interactions [2]. 

I start by examining thermodynamic fluctuations with a single variable, and em- 
phasize the problems that result on attempting to go beyond the Gaussian approx- 
imation. A corrected theory is given, and when a second independent fluctuating 
variable is added to this theory the meaning of R emerges. \R\ gives the correlation 
length and the sign of R corresponds to whether the interparticle interactions are 
effectively attractive or repulsive. 

2 ONE FLUCTUATING VARIABLE 

Consider a thermodynamic system consisting of a set of microscopic constituents 
contained in a box of volume V with an internal energy U. Imagine that we may 
specify the thermodynamics completely by writing the system's entropy S* as a 
function of U and V, 

S = SiU,V), (1) 

the fundamental equation pj. Examples for which S{U,V) is known exactly are 
a gas of electromagnetic radiation (photons) and a paramagnet of noninteracting 
spins. Although the fundamental equation is rigorous only in the thermodynamic 
limit of infinite volume, thermodynamic applications use Eq. ([I]) for systems of any 
volume. Because 5", U, and V are all additive, the fundamental equation may be 
written as 

S = Vs{u), (2) 

where {s,u} = {S/V, U/V} are quantities per volume [3]. 

Let this thermodynamics constitute a large, closed universe Aq with constant 
internal energy Uq and constant volume Vq — )■ oo. Imagine Aq to be separated into 
two parts joined by an immovable, open partition enclosing a finite system A, with 
fluctuating internal energy U and constant volume V. The rest of the universe is 
the reservoir Ar, with fluctuating internal energy Ur and constant volume K- — )■ cxo 
as illustrated in Fig. 1. Quantities with subscripts and r refer to properties of 
Aq and A^, respectively. From Eq. ^ we know that a single entropy per volume 
function s = s{u) is common to all three systems Aq, A^., and A. 



A and Ar exchange energy through fluctuations, which are described by thermo- 
dynamic fluctuation theory. The microcanonical ensemble asserts that all accessible 
microstates of Aq occur with equal probability |3]. Therefore, the probability of 
finding the internal energy per volume of A between u and u + du is proportional 
to the number of microstates of Aq corresponding to this range: 

P{u, V)du = Cno{u, V)du, (3) 

where ^o{u, V) is the density of states, and C is a normalization factor. Boltzmann's 
expression for the entropy S = fc^lnfi, where ks is Boltzmann's constant, yields 
Einstein's famous relation [^: 



P{u, V)du = Cexp 



Soiu,V) 



du, (4) 



where So{u,V) is the entropy of Aq when A has internal energy per volume u. 
Equation Q is the basis of thermodynamic fluctuation theory. 
Because entropy is additive, we have 

So{u,V)=Vs{u) + VrS{Ur). (5) 

In statistical mechanics we assume that there is some unique division of energy 
between A and Ar that maximizes So{u, V), corresponding to equilibrium. Let this 
maximum correspond to u = u* and Ur = u*, and define the differences Au = u — u* 
and Auj. = Ur—u*. We expand both entropy densities in Eq. ([s]) about the maximum 
and obtain 

So(m, V) = §0 + Vs\u*)Au + Vrs'{u;)Aur + V^^s"{u*){Auy 

+ Vry'{u:.){Aurf + Vy"{u*){Auf + ■ ■ • , (6) 

where the prime indicates differentiation with respect to u, and 5*0 is the entropy of 
Aq in equilibrium. 

Expansion about a maximum means that the first-order terms in Eq. ^ must 
add to zero. By conservation of energy, 

VAu = -VrAUr, (7) 

and the cancellation of the first-order terms requires that s'{u*) = s'{u*). This 
cancellation forces u* = u* because the inverse temperature s'{u) is assumed to 
decrease monotonically with increasing u. Additivity of energy further requires 



that Vu* + VrU* = VqUq, leading to u* = m* = uq- The differences now become 
{Au, Aur} = {u — Uo,Ur — Uo}. Energy conservation expressed in Eq. ([T]) also means 
that terms beyond first order in Eq. ^ from expanding s{ur) will be negligible 
compared with the corresponding terms from expanding s{u). Hence, 



So{u, V) = So + Vy'{uo){Auy + Vy"{u,){Auf + 



(8) 



Truncating this series at second order leads to the Gaussian approximation in 
thermodynamic fluctuation theory [H |5] : 



Pg{u, V)du 




exp 



~^g{uo){Auf 



g{uo) du, 



where 



g{uo) = 



-s"iuo) 



(9) 



(10) 



is always positive because it is proportional to the heat capacity. The normalization 
factor in Eq. ^ results from Gaussian integration, which also yields the averages 



(Am) = f Au Pg{u, V)du = 0, 



(11) 



and 



(Au) 



{Auf Pg{u, V)du 



Vgiuo) ' 



(12) 



Let's move on to some points not found in books on statistical mechanics. First, 
although the conservation rule (Au) = results almost trivially from the Gaussian 
integration Eq. (11), it is easy to construct an argument that it must hold for 
any physically correct probability density P{u,V) for all V [6j. However, beyond 
the Gaussian approximation, thermodynamic fluctuation theory does not obey the 
conservation rule. For decreasing V, fluctuations in Am increase by Eq. (12), and 
discarding the third-order term in the entropy expansion in Eq. (Isj) becomes in- 



creasingly difficult to justify. If in Eq. (11) we were to calculate the integral with 



P{u, V) containing a third-order term, we would find (Am) j^ 0, because the proba- 
bility density is no longer symmetric about u = Uq. (Am) 7^ 0, violates the required 
conservation rule. 

Another problem associated with going beyond the Gaussian approximation be- 
comes evident on changing the thermodynamic coordinate m to a general coordinate 



x{u). In the microstate counting argument which led to Eq. Q, no special proper- 
ties of u were used (conservation of energy conies only later in the argument). We 
could as well have used any other thermodynamic function x = x{u) and written 



P{x, V)dx = C*exp 



So{x,V) 



dx. 



(13) 



where C is a normalization factor. But, innocuous as it seems, Eq. (13) is inconsis- 
tent with Eq. ^. If we transform coordinate u to 
(|4]), we obtain 



x[u 



on the left-hand side of Eq. 



P{u, V)du 



P{u,V) 




dx = P{x, V)dx. 



(14) 



On the right-hand side of Eq. ^ the function Sq{u, V) = So{x, V) remains the same 
under this transformation because entropy is a function of state, but du = {du/dx)dx 
picks up the Jacobian J = {du/dx). Because J varies with the thermodynamic state, 
it cannot simply be absorbed into the normalization factor, and we end up with 



P{x, V)dx 




exp 



S,{x,V) 



k 



B 



dx. 



(15) 



which is inconsistent with Eq. (13). 



We call an inconsistency resulting from a change in coordinates a lack of co- 
variance [2]. It means that either Eq. dll) or Eq. (13) or both are incorrect. It is 
significant that this problem does not arise in the Gaussian approximation. If we 
change the coordinate in Eq. (pi), we obtain the transformed Gaussian approxima- 
tion for the X coordinate 



Pg{x, V)dx 




exp 



~g{xo){Axf 



g{xo) dx, 



(16) 



where the transformed function g{xo) is 



9{xo) = g{uo) 



' du 
dx 



(17) 



with Ax = X — xq and xq = x{uo). The Gaussian approximation is clearly covariant 
as long as we use the transformation rule Eq. (17). Of course, we require small 
fluctuations for its validity. 



Coming back to the u coordinate, observe that if g{u) were a constant function, 
then all derivatives of s{u) from third order and up would be zero, and the Gaussian 
approximation, which is covariant and obeys the conservation rule, is exact for all V. 
The function g{u) is given by the thermodynamics of the system of interest and thus 
we have no control over whether or not it is constant. But could we transform to 
another thermodynamic coordinate x = x{u) where the function g in the Gaussian 
approximation is constant, and higher order terms do not appear? Addressing this 
question is a major theme of this paper. But the answer is not to be found in Eq. 
(|4]), which is fundamentally unsound beyond the Gaussian approximation. We need 
a new approach. 

We start by adding a second, larger partition concentric with the first. This 
partition encloses a system Ai of which A is part as shown in Fig. 2. The infinite 
reservoir A,, is now everything in Aq outside of Ai. Because Ai is larger than A, 



fluctuations in Aui will be less than those of Am, by Eq. (12), and the Gaussian 
approximation should work better for it. We now calculate the fluctuation proba- 
bility of u about ui by imagining Ai to be in a state with energy per volume ui. 
Consider the reservoir surrounding A to be system Ai — A, rather than the infinite 
Aj.. Expand the entropy of Ai about its maximum as in Eq. ([6]). Again, the first- 
order terms drop out, but, in contrast to the case with an infinite reservoir, we must 
keep both quadratic terms. I omit the details here (see Ref. [2]). The key is that 
the resulting fluctuations in Au about ui are smaller than the corresponding ones 
about Mo, with the state of Ai unspecified, and the Gaussian approximation should 
thus be improved. 

Two probability distributions {A fluctuating about Ai and Ai fluctuating about 
Ao) are usually more than we ultimately need. Multiplying the two probability 
distributions and integrating out the intermediate state ui of Ai, gets us what we 
want: the probability of A fluctuating about A^. The result should be better than 
that resulting from Eq. ^ alone, because the pair of Gaussian approximations are 
each improved, and because in the integral with respect to Ui we can evaluate g{u) 
at the variable Ui rather than always at the fixed Uq. 

We can improve on this two probability structure by going to a continuum limit 
with more concentric partitions, filling Aq densely with them from very large down 
to A. Each system fluctuates about the state of the one immediately larger than it 
according to a Gaussian approximation. To obtain the probability of finding A in 
some state, given the state of Aq, we integrate out all the intermediate states. Such 
a multiple integration is called a path integral. To make it covariant and have it 
satisfy the conservation rule for all V requires great care. 1 skip the difficult details 
here [7], and approach this problem instead from another direction. 



Let's work on finding a partial differential equation for P{u, V). In many cases 
in physics, fundamental principles obey partial differential equations, and there is 
no reason why this would not be the case here. Finding an acceptable partial 
differential equation much limits the possible ways the modified thermodynamic 
fluctuation theory could be expressed. To this end, note first that the form of the 
Gaussian approximation Eq. ^ looks like the solution of a diffusion equation, with 
the role of "time" played by the inverse volume 

t^l^, (18) 

and with a Dirac delta function "initial condition" (t = 0) peaked at An = 0. I 
now explore this connection in a general coordinate x = x{u), and will return to the 
coordinate u later. 

The most general form of the diffusion equation which preserves normalization 
[8] may be written as [21 [9] 

1 d^ 



dP 
~dt 



d_ 
dx 



[K{x)P] 



2 9x2 



-1, 



x)P 



(19) 



where the drift term K{x) and g{x) are both functions of x to be determined. We 
understand that nothing in our problem is diffusing in time; t is simply a measure 
of volume, and the "diffusion" reflects the fact that the thermodynamic state of 
A becomes increasingly uncertain as V decreases, and the probability distribution 
for u gets broader, by Eq. (12). If we can determine K{x) and g{x), then they. 



along with the initial condition and appropriate boundary conditions, complete the 
covariant theory Pj. 

For constant functions K{x) = K and g{x) = g it is straightforward to verify 



that the exact normalized solution to Eq. (19) is the Gaussian expression 



PG{x,t) 



9 

2nt 



exp 



1 
'2t 



g{Ax - Ktf 



(20) 



where Ax = x — Xq, and Xq is a constant. This solution has the Dirac delta function 
S{x — Xo) start at t = 0. We can also readily verify the Gaussian integrals 



(Ax) = f AxPG{x,t)dx = Kt, 



and 



(Ax) 



{Axf PG{x,t)dx 



t{l+gKH) 



(21) 



(22) 



both valid for all t. 



For small t the Gaussian expression in Eq. (20) peaks sharply about the start 



ing value Xq, and PG{x,t) does not sample the functions K{x) and g{x) over any 
substantial range that deviates from xq. Taking K{x) and g{x) to be constant (and 
evaluated at xq) in Eq. (19) is thus justified for small t. Deviations of the true 



P{x,t) from the Gaussian expression is related to P{x,t) being spread over a large 
enough range that we can no longer consider K{x) and g{x) to be constant in solving 
Eq. 



(19) 



The Gaussian integrals above connect both K and g to the thermodynamics 
of the system under consideration. Take t to be small, and switch back to the u 



coordinate. Compare Eqs. (20), (21), and (22) with the corresponding ones from 
the thermodynamic Gaussian expression Eqs. ([9]), (11), and (12). Observe that 
we must have K = and t = 1/V. Also observe that g is common to the two 
cases, is given by Eq. (10), and determines the width of the fluctuations. At flrst 



sight it might appear that the function K(x) in Eq. (19) is irrelevant because in u 
coordinates K{u) = 0. But nonzero K{x) turns out to be necessary to transform 
coordinates correctly. 

To make further progress with the theory, we must learn how to transform Eq. 
(19) to another coordinate x = x{x). First we must ask, is it even possible to 
transform P{x,t), K{x 
equation is unchanged: 

dt ~ dx \- ^''' \ ^ 2 dx^ 



and g{x) such that the form of the partial differential 



g-U. 



{x)P 



(23) 



That the partial differential equation keeps its form is required for covariance. Be- 
cause probability is a scalar quantity independent of the choice of coordinate, we 
require Pdx = Pdx, and thus 



We substitute this expression into Eq. (23), replace d/dx by {dx/dx)d/dx, divide 



by dx/dx, and equate coefficients of the corresponding derivatives of P with those 



in Eq. (19). The result is equality if and only if the functions K and g transform as 

f dx\ 
dx I 



9 = 9 



(25) 



and 



K 



(JjUj -L -f (Jj Uj 

dx 2 dx'^ 



(26) 



Hence, the partial differential equation Eq. ( 19 ) is covariant if we use the appropriate 
transformation rules. I add that it was not a priori clear that we could find acceptable 
transformation rules allowing covariance. Our success further attests to the validity 



of Eq. (19). The transformation rules and the connection to thermodynamics in the 
u coordinate determine K{x) and g{x) uniquely for any coordinate x = x{u). We 
have thus solved the problem of fluctuation theory with one variable. 

Is it possible to transform to a coordinate x = x{u) with transformed K{x) and 



g{x) constant everywhere? If so, the Gaussian expression Eq. (20) would apply for 



all t. First observe that finding x = x{u) with constant g is always possible. By Eq. 



(25) we see that 



^ = y9 V V^(") ^"' (27) 



with any constant g, does the job. However, knowing x{u) leaves us with no more 
degrees of freedom, and we generally cannot also solve Eq. ( |26[ ) for constant K. Our 
attempt at a perfect coordinate is only half successful. But notice that for small t it 
is g which determines the width of the fluctuations; K is needed only as a correction 
term to keep the average energy fixed as t becomes bigger (and is not even relevant 
for small t). Thus, although we have solved only half the problem, we have solved 
the better half. Regardless, a necessary condition for the Gaussian expression to 
hold for all t, constant g, may always be satisfied with one variable. 

I mention in passing another theme in addition to covariance and conservation 



of energy. There is more than one way to obtain P{x,t) from Eq. (19), and the 
various ways must yield the same results - the requirement of consistency. The first 
(and best) way is to solve for P{x, t) in a single step, starting from a delta function 
initial condition at t = 0. Also possible is a two step solution method, as in Fig. 
2, where we first solve for P{xi,ti) from t = to some ti < t for the probability 
of finding xi. Then we solve to t for the conditional probability P{x,t\xi,ti), the 
probability that we get x at t given that we know it is xi at ti. Integrating out the 
intermediate variable yields P{x,t) = J P{x,t\xi,ti)P{xi,ti)dxi. To be consistent, 
P{x,t) must agree with that from the single step method. 



I point out that solutions to Eq. ( 19 ) satisfy covariance, conservation, and consis- 
tency [2j, in contrast to P{x,t) calculated from classical thermodynamic fluctuation 
theory, Eq. (Ill), which satisfies none of these conditions beyond the Gaussian ap- 



proximation. I claim that the theory based on Eq. (19) is the correct way to extend 
thermodynamic fluctuation theory beyond the Gaussian approximation fUJl [TTl [T2] . 
Ruppeiner [2] has given an explicit calculation for the simple paramagnet, including 
the issue of boundary conditions (not considered here). The covariant theory shows 



considerable improvement over the classical one. 

The improvements resulting from the covariant theory with just one variable 
affect little in physics, because problems characterized by fundamental equations 
of the form S = S{U, V) are generally well understood. But, the story changes if 
we add another variable. A fundamentally new object appears, the thermodynamic 
curvature. 



3 TWO FLUCTUATING VARIABLES 

Consider a pure fluid system (gas or liquid) of A^ identical particles in a volume V 
which may be open or closed. The fundamental equation is |3j 

S = S{U,N,V). (28) 

Because S, U, N, and V are all additive, we may write instead 

S = Vs{u,p), (29) 

where {s, u, p} = {S/V, U/V, N/V} are quantities per volume. 

We again partition a closed, infinite universe Aq into two parts, a finite A and 
an infinite reservoir A^, by a fixed, imaginary partition as shown in Fig. 1. Both U 
and N, but not V, fluctuate |3]. Define {a^jO^} = {u,p}, and 

F. . ^. (30) 

where a = 1,2. Basic thermodynamics gives {^1,^2} = {l/T,—p/T}, where T is 
the temperature, and p is the chemical potential. The properties of Aq and Ar are 
denoted by subscripts and r, respectively. 

The probability of finding the state of A in the range (u, p) to {u + du, p + dp) is 



Pdudp = Cexp ( - — ) dudp, (31) 

KKrJ 



' Sq 
where C is a normalization factor, and 5*0 may be written as 

5*0 = Vs{u,p) + VrS{Ur.,Pr)- (32) 

If we apply the same argument as used in Sec. II, there will be a maximum for Sq 



in the state where a"^ = a'^ = a^. I expand the entropies in Eq. (32) about this 
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maximum in powers of the differences Aa" = a" — Oq and Aa" = a" — Oq and obtain 

1 OF 

ASo = VF.Aa^ + VrFr,Aa^ + V-^Aa^^Aa'' 

where A 5*0 is the difference between 5*0 and its maximum value, and summation 
over repeated indices is assumed. 

Conservation of both energy and particles requires VAa"" = — K-Aa". A neces- 
sary condition for maximum entropy is thus Fa = Fra- Also, the second quadratic 



term in Eq. (33) is negligible compared with the first, and Eq. (33) becomes 



^ = -^^^,Aa'^Aa^ (34) 

Kb I 

where the symmetric intensive matrix 

must be positive-definite because ASq < 0. (This property holds only for short-range 
interactions typical of most fluid and magnetic systems. If, for example, gravity is 
present, it might be difficult to find a local entropy maximum.) 



If we substitute Eq. (34) into Eq. (^31| and calculate the normalization factor, 



we obtain the Gaussian approximation |1] 

Pcda^da^ = ( — j exp ( g^yAa^Aa'^j^da^da^^ (36) 

where g is the determinant of the matrix gap. 

Let's introduce a new component into our discussion and note that the positive- 
definite quantity 

Ae = g^^Aa^Aa" (37) 

has the appearance of a line element or distance between neighboring thermody- 
namic states [L3\, or more technically, a Riemannian metric, as described in the 
Appendix. In thermodynamics, distance clearly measures probability: the less the 
probability of a fluctuation between two states, the further apart they are. 
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To transform to another pair of thermodynamic coordinates (x^,x^), we note 
that because entropy is a function of state, A£^ keeps its value regardless of what 
coordinates in which it is we expressed. Hence, we can write 

Pg dx^dx^ = (— ) exp ( g^^Axf^Ax") ^ dx^dx^, (38) 

provided that we transform the metric elements gai3 as a second-rank tensor, as in 



Eq. (65). The tensor transformation property may look novel in thermodynamics. 



but it is at least implicit in Ref. [Ij. In general coordinates, the line element is 

Af = Qf^^Ax^Ax". (39) 

The Gaussian approximation yields the averages 

(Ax") = Iax'^Pg dx^dx^ = 0, (40) 

and 

(Ax^Ax^) = [Ax^Ax^Pg dx^dx"^ = ^, (41) 

where g""^ denotes the components of the inverse of the matrix gai3- 

Particularly simple are coordinates that diagonalize gap- For example, in (T, p) 
coordinates [1]: 

This line element is readily expressed in terms of the Helmholtz free energy per 
volume 

(43) 



(44) 



/ = 


u — Ts. 


We have 




s = 


(df\ 


and 




/^ = 


(df\ 

\9pJt' 



(45) 
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Equation (41) yields, 



(Ap)' 



V 




For the ideal gas, 

/ = pksT In p + pkBh{T), 



(46) 



(47) 



where h{T) is some function of the temperature with negative second derivative to 
assure a positive heat capacity. Equation (46) yields ((AA^)^) = (N). Integrating 



AT out in Eq. (38) gives the Gaussian approximation for fluctuations in AA^. This 



approximation may be compared with the exact Poisson distribution for A^. Results 
for (A^) = 5 are shown in Fig. 3 [iQl . For the ideal gas, the Gaussian approximation 
is effective even at size scales containing just a few particles. 

The addition of a second independent fluctuating variable again raises issues of 
covariance, conservation, and consistency. As in Sec. II, it is easy to construct the 
correct theory around an appropriate partial differential equation j2]: 



dP 

~dt 



d 1 d^ 



(4^ 



In {u, p) coordinates, K"^ 
formations are 



and gas is given in Eq. (35). The coordinate trans 



P = P 



dx 



dx 



(49) 



where \dx/dx\ is the Jacobian of x{x), 



dx" dx^ 



Qt^v - 9al3^^ 



dx^^ dx'' 



and 



K' 



dx" 
dx^^ 



K^ + L 



2" dxt'dx"' 



(50) 



(51) 



For constant K" and constant metric elements 9^/3 we can write down the Gaus- 
sian expression 

1 



PG{x,t) 



(2^) ^^"P 



2t 



^^^(Ax^ - K^t){Ax'' - K^t) 



(52) 
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which satisfies Eq. (48) exactly. Clearly, for small t this Gaussian expression will 
always be adequate since K°' and ga/s don't get a chance to fiuctuate very much. 
For larger t, with its larger fiuctuations, the key to getting the Gaussian expression 
to work is finding a transformation to coordinates with constant if" and Qap- As in 
Sec. II, a necessary requirement for this coordinate transformation is that it yield 
constant Qap- However, as I will argue in Sec. V, such a coordinate transformation 
does not generally exist, in contrast to the case with one fiuctuating independent 
variable. There will thus be a lower limit on V where the Gaussian expression must 
fail regardless of choice of coordinates. 

To conclude this section, I point out that Weinhold [1^ originated thermody- 
namic metrics in the form of inner products based on the Hessian of the internal 
energy. The positive-definite nature of this inner product represents the second 
law of thermodynamics. But Weinhold's geometry lacks a true Riemannian metric 
structure since it had no underlying physical notion of distance, such as is offered 



by fiuctuations in Eq. (37) where the Hessian of the entropy come into play [T3] . 
An entropy metric was also used by Andresen, Salamon, and Berry [15J as a mea- 
sure of the dissipated availability in finite-time thermodynamics. Ruppeiner |TOl Hi] 
first suggested using a covariant and consistent diffusion type equation for thermo- 
dynamic fluctuation theory. Diosi and Lukacs [12j completed this theory, adding 
conservation and drift, and were the first to write the equation in the form in Eq. 



(48). 



4 A PHYSICAL LIMIT 

As an application, consider a pure fiuid possibly near its critical point. The attractive 
part of the force between particles causes particles to be closer to each other on 
the average than they would be in a purely random situation. A clustered pair 
will attract a third particle, and so on, resulting in a droplet of particles of some 
characteristic size, the correlation length ^ j5|. This length is zero for the ideal gas 
and approaches infinity at the critical point. Unless we are in a regime where the 
fiuid coexists in two phases (assumed not to be the case), these droplets continually 
disassociate and reform. 

This picture is an oversimplification. The density of the fiuid is a function of 
position p = p{r), which has contour surfaces of constant density. Consider a surface 
of constant po- According to Widom |16], this surface is enormously complex and 
sponge-like. Any straight line makes many intersections with it, as shown in Fig. 4 
(reproduced from Widom [IS]). The mean distance between these intersections is a 
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measure of ^. Because these intersections occur independent of the direction of the 
hne, it is convenient to think of the regions p < po and p > po as each consisting 
of many small, independent volume elements ("droplets") of linear dimension ^ (see 
Fig. 4) though of course these volume elements are in reality connected in one 
sponge-like mass. 

The correlation length ^ is commonly assumed to require statistical mechanics 
for its proper definition and calculation. The key argument of this paper is that ^ 
can be deduced from thermodynamics. Notice in Fig. 4 that if ^ ^ ^^, then system 
A would see itself mostly surrounded by a local environment with density near to 
Po at any time. Being able to regard the reservoir Ar of A as a uniform system 
with density po is an essential requirement in the Gaussian expression, and it is thus 
justified for V ^ C,^. If, however, V <^ C,^, A would at any time most likely be in a 
region between the curves of constant density po in Fig- 4. The local surroundings 
now seen by A at any time has a thermodynamic state different from that of the 
overall system Aq. The density fluctuations would be typically bimodal, as in Fig. 
5. 

The length scale ^ thus marks the intrinsic breakdown of the Gaussian expression. 
In Sec. [5| I relate this breakdown to the mathematical structure we have been 
developing. 

5 THERMODYNAMIC CURVATURE 

For V sufficiently large, the physical limit discussed in Sec. |4]does not constrain us, 
and in {u,p) coordinates the Gaussian expression should work just fine. For small 
fiuctuations, the matrix elements gai3 will be nearly constant over the full range of 
fiuctuations. 

But even in this apparently safe volume regime we can always invalidate the 
Gaussian expression by transforming to "bad" coordinates in which g^i^ varies sharply 
for even small fiuctuations. However, this mathematical game has no relation to any- 
thing physically fundamental. We should simply avoid such coordinates. It is far 
more interesting to transform to "good" coordinates with the Gaussian expression 
extended to small V. As we have seen, a necessary property of good coordinates is 
slowly varying ga/s as A fiuctuates. Better yet, might we find "great" coordinates 
with the Gaussian expression extended to near atomic volumes? A reason to hope 
for "great" are the ideal gas results shown in Fig. 3, where we lucked into such 
coordinates. But the ideal gas, with ^ = 0, does not suffer from the physical break- 
down of Sec. |4l and thus it is no surprise that the Gaussian expression works to 
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such small volumes. More generally, ^ marks the physical onset of non-Gaussian 
fluctuations in all coordinates. I will now argue that this volume can be deduced 
from thermodynamics. 

The key is the Riemannian curvature scalar R of the metric Qaf^- Imagine V 



large enough for the Gaussian expression Eq. (52) to be valid in some choice of 
coordinates. For the Gaussian expression to be valid, we require constant gap over 
typical fluctuations. With constant gap, simple hnear algebra allows a transforma- 



tion to a locally fiat Cartesian line element, as in Eq. (69). Clearly, a valid Gaussian 
expression requires the system to be in the locally fiat geometric regime. 

This relation is elementary. Much more interesting are larger fluctuations. Define 
the squared distance 



r 



2 



g^^Ax^Ax^ (53) 



to be the size of a typical fluctuation of the state of A about that of A^. The 



argument leading to Eq. ( 70 ) leads us to the observation that 



r' - ^, (54) 

\ri\ 

beyond which we are no longer in the geometric locally flat regime. Beyond this 
approximate value of r, no local coordinate system with constant gap exists, and the 
Gaussian expression fails in all coordinate systems. 



The argument of the exponential in Eq. (38) is roughly —1 for typical fluctua- 
tions, corresponding to 

r^ ~ |,. (55) 



If we combine Eq. (55) with Eq. (54), we see that the Gaussian expression must 



fail in all coordinate systems if 

V <~ ^. (56) 

6 

This breakdown is intrinsic to the geometry, and has nothing to do with pick- 
ing bad coordinates. Nor may it be solved by any good coordinates. It calls for 
some physical interpretation. I relate it to the physical breakdown of the Gaussian 
expression in Sec. IV and interpret \R\ as being proportional to the correlation 
volume 

\R\ ~ e", (57) 

where d is the dimensionality of space. 
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6 CALCULATION OF R 

Given some statistical mechanical model, we may test the proposition that \R\ ~ ^'^ 
by calculating both R and ^. We first evaluate R for the ideal gas and then follow 
with a brief review of numerous other models for which R has been evaluated. I also 
make some observations about the sign of R. 

Although R for the ideal classical gas requires almost no calculation, this cal- 
culation nevertheless makes a good exercise. We already know \R\ must be small 
because if particles move randomly with respect to each other, there is no tendency 



for them to cluster. Equations (42)- (45) and (47) yield the line element in coordi- 
nates {x^,x'^} = {T,p}: 

Ae^ = -P!^AT^ + -Ap^ (58) 

T P 

from which we can read off the metric elements 



We substitute these metric elements into Eq. (68) for R and find the exact result 



R = (60) 

for the ideal classical gas jT3] . 

This result came as a surprise, and pointed to some mysterious connection be- 
tween curvature and interactions. But the argument given in Sec. V clears up 
the mystery. Also initially surprising but now understood is the proportionality 



\R\ ~ ^'^ in Eq. (57). It has been confirmed for several statistical mechanics models 
[ISlIIZlIIHlIIglEOlEIlEaEalEllES], particularly near critical points, where ^'^ is 
large enough to encompass a volume containing many particles. I add that R is typ- 
ically negative in my sign convention (see Appendix) in the standard critical point 
models, for example, van der Waals for which the attractive interaction dominates 
the long-range behavior. 

Another key example is an ideal gas of bosons. This system has no interparticle 
interactions, but quantum statistics leads particles to bunch closer than in a classical 
ideal gas. This clustering should be measured by R. The ideal Bose gas suffers a 
Bose-Einstein condensation where a macroscopic number of particles pool in the 
ground state at low temperature. It was found that R is negative and diverges to 
negative infinity at absolute zero [26] . 
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Another example is an ideal gas of ferniions. R was found to be positive and 
to diverge to positive infinity at absolute zero [261 [27]. Janyszek and Mrugala [26] 
have emphasized the difference in the sign of R between the Bose and Fermi gases. 
An explicit connection with the correlation length for ideal quantum gasses remains 
to be estabhshed. 

Table [I] shows the sign and divergence of R for several systems [28]. R is negative 
where attractive interactions dominate, and positive R where repulsive interactions 
dominate [22] • Table [I] includes several weakly interacting systems with R on the 
order of the volume of an interparticle spacing or less. I interpret such values of R 
as physically equivalent to zero, because the meaning of ^'^ of this size is lost in the 
noise as thermodynamics breaks down as individual particles or spins become visible. 
The one-dimensional antiferromagnetic Ising model [TTJ [31] may be misplaced in this 
category, since the ordering field in antiferromagnetic systems is a staggered field, 
and not the constant field used for the calculations in Table [H 

There are four systems in Table [l] for which R is both positive and negative. 
Not all these have been analyzed in terms of attractive and repulsive interactions. 
The sign of R for the one-dimensional g-state Potts model [211 125] is related to 
q. For q > 2 and nonzero magnetic field, there are significant regions of negative 
R at low temperatures. The Takahashi gas [TS] has negative R in the gas-like 
phase where attractive interactions dominate, and small |-R| in the liquid-like phase 
where interactions are effectively short-range. Going from one phase to the other 
by changing the density at constant temperature gives a pseudophase transition 
accompanied by a sharp positive spike in R. An abrupt change in the sign of R is 
also present in the one- dimensional Ising ferromagnet of finite N spins [21]. R is 
negative for large A^, and sharply increases to large positive values as A^ is decreased. 
Both positive and negative signs for R are present in the noninteracting gas with 
Gentile's statistics [2Z]. Such a gas allows a maximum number of particles p into any 
quantum state. In the classical limit \R\ is always small, R is negative for p > 2, and 
R is positive for the Fermi gas {p = 1). In the degenerate limit T — )■ 0, i? — )■ +oo 
for all finite values of p. Oshima et al. [2Z] also corrected some earlier Fermi gas 
results [26] . Also in the category of having R with both signs is a gas of anyons [35] , 
particles in two dimensions with fractional spin. This model has a parameter whose 
variation allows us to change it continuously from a Bose gas to a Fermi gas; the 
sign of R changes correspondingly from negative to positive. 



7 CONCLUSIONS 

Can thermodynamic curvature tell us anything we could not obtain otherwise? If we 
have a microscopic model, statistical mechanics offers powerful methods for calculat- 
ing its thermodynamic properties. It also usually allows us to calculate microscopic 
properties, such as the correlation length. Calculating R might not teach us much 
about such a model. 

If the system's underlying microscopic structure is unknown, but we know its 
thermodynamic properties, calculating R gives a bigger payoff. Such a case occurs 
for black holes which has a well developed thermodynamics [37]. Little is known 
about its possible microscopic foundations. There have been several calculations of 
the thermodynamic curvature R for black holes [281 EH] ■ It has been found [28] that 
the sign of R is almost always positive for the spinning, charged Kerr-Newman black 
hole, suggesting that an appropriate microscopic model might consist of some type 
of repulsive fermions rather than one of the standard critical point models listed in 
Table I. 

If R does augment thermodynamics in some fundamental way, then we must 
investigate why and to what extent. For example, we have interpreted just \R\. 
Beyond the observations summarized in Table [1} it is not clear what the sign of R 
signifies and why. 

\R\ relates to hyperscaling in critical phenomena [5], [16]. Hyperscaling relates an 
appropriate free energy to the correlation length by </> ~ ^~'^. Combining this with 
\R\ ~ ^'^ yields \R\ ~ (j)~^ [SU]- Because R may be expressed in terms of (p and its 
derivatives, we have a partial differential equation for 0. Its solution, in conjunction 
with boundary conditions, yields using no statistical mechanics at all. Tests have 
been few, so any review of results here would be premature. 

8 Appendix: Riemannian geometry 

Consider a two-dimensional surface embedded in a three-dimensional Euclidean 
space with Cartesian coordinates {x,y,z) via some function z = z{x,y). For ex- 
ample, a two-sphere centered at the origin with radius a gives z = iA/a^ — x^ — y'^, 
as in Fig. 6. 

The two-dimensional surface, with some choice of coordinates (x^,a;^), inherits 
its distance from the surrounding Euclidean space. Its line element is [10] 

Af = g^^Ax^'Ax", (61) 
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where /z, z/ = 1, 2, and the metric elements gai3 constitute a positive-definite matrix. 
For example, the two-sphere in spherical coordinates {6, 0) has 

Af = a^Ae^ + a" sin^ 9 A(p\ (62) 

The distance between pairs of neighboring points does not depend on the choice 
of coordinate system, and hence A(? transforms as a scalar. Consider some other 
pair of coordinates (x^,x^) in which the metric elements are cjap- Because 

Ax- = ^Ax'^, (63) 

we have 

^^' = ^-'rn^W^^^'^^'^ " ^M.A.x'^AS^ (64) 

The transformed metric elements are thus 
dx^ dx 



Equation ( 65 ) is the transformation rule for a second-rank tensor [iQl HI] . The two- 
dimensional surface also has area element y/gdx^dx"^ which transforms to y/gdx^dx'^, 
with g and g the determinants of gaj3 and gap, respectively [30] . 

We measure distance along any curve on the two-dimensional surface by inte- 
grating the line element along it: 



g^^dxi^dx'^. (66) 

Introduce curves on the two-dimensional surface minimizing the distance between 
pairs of points. These curves are called geodesies. Imagine a simple exercise starting 
from some central point anywhere on the two-dimensional surface and drawing along 
geodesies a fixed distance r (not too large) in all directions. The resulting endpoints 
define a circle; see Fig. 6 for the circle on the two-sphere. 

Locally, any two-dimensional surface looks fiat, and a very small circle has cir- 
cumference C ~ 27rr. Draw out further and deviations from locally Euclidean 
geometry appear. For example, for the circle on the two-sphere in Fig. 6, C is less 
than 27rr because the circle lies entirely on one side of a plane tangent to the two- 
sphere at the center of the circle. A challenging exercise in non-Euclidean geometry 
shows that 



C = 27rr + -Rr^ + 0(r'^), (67) 

6 
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where the Riemannian curvature scalar R gives the size of the leading correction to 
Euclidean geometry. 

The value of R depends on where we are on the surface. (Historically, the 
Gaussian curvature K = —R/2 came first. K is unambiguously defined and is 
positive for the two-sphere. The sign of R is subject to convention. I use Weinberg's 
sign convention [30j, with R negative for the two-sphere.) Because R is based on 
lengths of curves, it is independent of the choice of coordinate system, and is a 
purely geometrical quantity. For a plane R = 0. For the two-sphere in Fig. 6 a 
simple calculation shows R = — 2/a^. For a diagonal metric we have 



\f9 




1 dg-n 




(68) 



Surfaces of constant positive R also exist, pseudospheres as in Fig. 7 |12]. Here, 
a small circle drawn out from any point lies on both sides of a plane tangent to the 
two-sphere at the center of the circle, and C is larger than that of a circle with the 
same radius drawn on a plane. 

For simplicity, I have emphasized drawing small circles. Other questions may be 
asked regarding deviations from Euclidean geometry. For example, how much does 
the sum of the interior angles of a triangle with geodesic sides deviate from vr? R 
answers all such questions. 

In a small neighborhood of a non-singular point, the geometry is always flat. 
Hence, we may find a pair of locally Cartesian coordinates (X, Y) on the two- 
dimensional surface with 

A£2 = AX^ + Arl (69) 

How large may r become before "fiat" becomes inadequate? For a rough estimate. 



set the correction term in Eq. (67) equal to the leading term, yielding 

r' - ^. (70) 

\it\ 

For r beyond this limit, the geometry is intrinsically non-Euclidean, and may not be 
represented with locally Cartesian coordinates because they will always give 27rr for 
the circle circumference. Neither is it possible to find any local coordinate system in 
which we can treat gap as effectively constant. If such a coordinate system existed, 
we could transform these coordinates to locally Cartesian coordinates, which we 
know cannot represent the curvature of the surface. 

The discussion has been of the surface geometry originated by Gauss [13] which 
starts with an embedding function z{x,y), and a surface metric induced by the 
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embedding in the three-dimensional Euchdean space. The approach in thermody- 
namics is the converse. We start with a surface metric, and seek a corresponding 
embedding function. But this converse problem usually has no solution, and an em- 
bedding function rarely exists [H]. Therefore, Gauss's surface geometry is in most 
cases useless in metric thermodynamics. 

To solve this problem, we need to go one step further to Riemann's geometry. 
Riemann kept as basic the notion of a smooth, two-dimensional surface of points with 
coordinates (a;^,x^). (More technically, a two-dimensional manifold. Generalization 
to higher dimensions is possible.) He dropped the embedding in a three-dimensional 



Euclidean space, and with it the function z{x, y), but kept the line element Eq. (61 ) 



lengths of curves Eq. (66), and deviations from Euclidean geometry characterized 



by measurements within the two-dimensional surface, such as those leading to Eq. 



(67). Riemann also kept equations independent of embedding, such as Eq. (68) 
for R. Thermodynamic metric geometry requires a switch to Riemann's picture. 
However, for visualizations and calculations we are usually safe using the ideas and 
many relations from Gauss' geometry. 

I thank Harvey Gould for many suggestions which improved the narrative flow. 



References 

[1] A. Einstein, "On the general molecular theory of heat," Ann. Phys. (Leipzig) 14, 354-362 
(1904); ibid., "The theory of the opalescence of homogeneous fluids and liquid mixtures near 
the critical state," 33, 1275-1298 (1910). 

[2] G. Ruppeiner, "Riemannian geometry in thermodynamic fluctuation theory," Rev. Mod. 
Phys. 67, 605-659 (1995); "Erratum: Riemannian geometry in thermodynamic fluctuation 
theory [Rev. Mod. Phys. 67, 605 (1995)]," Rev. Mod. Phys. 68, 313-313(E) (1996). 

[3] H. B. Callen, Thermodynamics and an Introduction to Thermo statistics (John Wiley & Sons, 
New York, 1985). A gas of electromagnetic radiation with S = S{U,V) is discussed in Sec. 
3-6. 

[4] L. D. Landau and E. M. Lifshitz, Statistical Physics (Pergamon, New York, 1977). 

[5] R. K. Pathria, Statistical Mechanics (Butterworth-Heinemann, Oxford, 1996). 

[6] For this argument, divide Aq up into 77 identically shaped disjoint systems Ai, A2, ..., A^ each 
of the same volume V, and joined by open partitions. At some time, system Ai contains energy 
Ui and has energy density Ui — Ui/V — rjUi/Va. Additivity of energy yields uq = T,Ui/r]. 
Averaging both sides, noting that Uq is constant, and taking (ui) the same for all i, establishes 
(Am) = 0. If we worry about systems near the boundary of Aq being special, impose periodic 



22 



boundary conditions on Aq. In the thermodynamic hmit Vb — > oo, the choice of boundary 
conditions doesn't matter. 

[7] H. Grabert and M. S. Green, "Fluctuations and nonhnear irreversible processes," Phys. Rev. 
A 19, 1747-1756 (1979). The full treatment is given here, with an arbitrary number of 
variables. The physics deals with irreversible thermodynamics, and so is diiferent from that 
in this paper. But the mathematical reasoning should be analogous. 



We may easily verify preservation of normalization by integrating both sides of Eq. ( 19 1 with 



respect to x and setting P{x, t) and its first derivative to zero at the limiting values of x. 

[9] R.Graham, "Covariant formulation of non-equilibrium statistical thermodynamics," Z. Phys. 
B 26, 397-405 (1977). 

[10] G. Ruppeiner, "New thermodynamic fluctuation theory using path integrals," Phys. Rev. A 
27, 1116-1133 (1983). 

[11] G. Ruppeiner, "Thermodynamic critical fluctuation theory?," Phys. Rev. Lett. 50, 287-290 
(1983). 

[12] L. Diosi and B. Lukacs, "Covariant evolution equation for the thermodynamic fluctuations," 
Phys. Rev. A 31, 3415-3418 (1985). 

[13] G. Ruppeiner, "Thermodynamics: A Riemannian geometric model," Phys. Rev. A 20, 1608- 
1613 (1979). 

[14] F. Weinhold, "Thermodynamics and geometry," Phys. Today 29 (3), 23-30 (1976). 

[15] B. Andresen, P. Salamon, and R. S. Berry, "Thermodynamics in finite time," Phys. Today 
37 (9), 62-70 (1984). 

[16] B. Widom, "The critical point and scaling theory," Physica 73, 107-118 (1974). 

[17] G. Ruppeiner, "Application of Riemannian geometry to the thermodynamics of a simple 
fluctuating magnetic system," Phys. Rev. A 24, 488-492 (1981). 

[18] G. Ruppeiner and J. Chance, "Thermodynamic curvature of a one-dimensional fluid," J. 
Chem. Phys. 92, 3700-3709 (1990). 

[19] D. Brody and N. Rivier, "Geometrical aspects of statistical mechanics," Phys. Rev. E 51, 
1006-1011 (1995). 

[20] B. P. Dolan, "Geometry and thermodynamic fluctuations of the Ising model on a Bethe 
lattice," Proc. R. Soc. Lond. A 454, 2655-2665 (1998). 

[21] B. P. Dolan, D. A. Johnston, and R. Kenna, "The information geometry of the one- 
dimensional Potts model," J. Phys. A: Math. Gen. 35, 9025-9035 (2002). 



23 



[22] W. Janke, D. A. Johnston, and R. P. K. C. Malmini, "Information geometry of the Ising 
model on planar random graphs," Phys. Rev. E 66, 056119-1-5 (2002). 

[23] W. Janke, D. A. Johnston, and R. Kenna, "Information geometry of the spherical model," 
Phys. Rev. E 67, 046106-1-4 (2003). 

[24] D. C. Brody and A. Ritz, "Information geometry of finite Ising models," J. Geom. Phys. 47, 
207-220 (2003). 

[25] D. A. Johnston, W. Janke, and R. Kenna, "Information geometry, one, two, three (and four)," 
Acta Phys. Pol. B 34, 4923-4937 (2003). 

[26] H. Janyszek and R. Mrugala, "Riemannian geometry and stability of ideal quantum gases," 
J. Phys. A: Math. Gen. 23, 467-476 (1990). 

[27] H. Oshima, T. Obata, and H. Hara, "Riemann scalar curvature of ideal quantum gases 
obeying Gentile's statistics," J. Phys. A: Math. Gen. 32, 6373-6383 (1999). 

[28] G. Ruppeiner, "Thermodynamic curvature and phase transitions in Kerr-Newman black 
holes," Phys. Rev. D 78, 024016-1-13 (2008). 

[29] There is no connection between the sign of R and thermodynamic stability. The thermody- 
namic metric elements ga/s always constitute a positive-definite matrix for ordinary fluid and 
magnetic systems, regardless the sign of R. The Bose and Fermi gasses both illustrate this 
property. A good, purely geometric, example is offered by the sphere and the pseudosphere 
in the Appendix. 

[30] S. Weinberg, Gravitation and Cosmology (John Wiley & Sons, New York, 1972). 

[31] H. Janyszek and R. Mrugala, "Riemannian geometry and the thermodynamics of model 
magnetic systems," Phys. Rev. A 39, 6515-6523 (1989). 

[32] D. C. Brody and D. W. Hook, "Information geometry in vapour- liquid equilibrium," J. Phys. 
A: Math. Thcor. 42, 023001-1-33 (2009). 

[33] G. Ruppeiner, "Equations of state of large gravitating gas clouds," Astrophys. J. 464, 547- 
555 (1996). 

[34] G. Ruppeiner and C. Davis, "Thermodynamic curvature of the multicomponent ideal gas," 
Phys. Rev. A 41, 2200-2202 (1990). 

[35] B. Mirza and H. Mohammadzadeh, "Nonperturbative thermodynamic geometry of anyon 
gas," Phys. Rev. E 80, 011132-1-9 (2009). 

[36] K. Kaviani and A. Dalafi-Rezaic, "Pauli paramagnetic gas in the framework of Riemannian 
geometry," Phys. Rev. E 60, 3520-3525 (1999). 

[37] J. D. Bekenstein, "Black-hole thermodynamics," Phys. Today 33 (1), 24-31 (1980). 

24 



[38] J. E. Aman, J. Bedford, D. Gmmillcr, N. Pidokrajt, and J. Ward, "Ruppeiner theory of black 
hole therraodynamics," J. Phys.: Conference Series 66, 012007-1-10 (2007). 

[39] G. Ruppeiner, "Riemannian geometric theory of critical phenomena," Phys. Rev. A 44, 
3583-3595 (1991). 

[40] G. B. Arfken and H. J. Weber, Mathematical Methods for Physicists (Academic Press, New 
York, 2001). 

[41] D. Laugwitz, Differential and Riemannian Geometry (Academic Press, New York, 1965). 

[42] C. Hsiung, A First Course in Differential Geometry (John Wiley & Sons, New York, 1981). 

[43] K. Gauss, General Investigations of Curved Surface (Raven, Hewlett, New York, 1965). 

[44] B. Andresen, R. S. Berry, R. Gilmore, E. Ihrig, and P. Salamon, "Thermodynamic geometry 
and the metrics of Weinhold and Gilmore," Phys. Rev. A 37, 845-848 (1988). 



25 



Table 



System 


n 


rf 


Sign of R 


Divergence 


3D Bose gas 26] 


2 


3 


— 


T^O 


ID Ising ferromagnet 1 7| 131] 


2 


1 


— 


T^O 


Critical region 21 [131 [IS] 


2 




— 


critical point 


Mean-field theory 31 


2 




— 


critical point 


vanderWaals 21 [lyl|32j 


2 


3 


— 


critical point 


Ising on Bethe lattice 20] 


2 




— 


critical point 


Isine on 2D random graph 22[ 125] 


2 


2 


— 


critical point 


Spherical model [231 [2S 


2 




— 


critical point 


Self-gravitating gas J33] 


2 


3 


— 


unclear 


ID Ising antiferromagnet |17l 131 


2 


1 


— 


R\ small 


Tonks gas \m\ 


2 


1 


— 


R\ small 


Ideal gas 13] 


2 


3 





R\ small 


Ideal paramaenet 17^ 131] 


2 







R\ small 


Multicomponent ideal gas ^l 


>2 


3 


+ 


R\ small 


ID Potts model 211 EH] 


2 


1 


± 


T^O 


Takahashi gas 1 S 


2 


1 


± 


T^O 


Finite ID Ising ferromagnet j24j 


2 


1 


± 


T^O 


Anyon gas 35 


2 


2 


± 


T^O 


2D Fermi gas 2E] 


2 


2 


+ 


T^O 


3D Fermi gas [26l [27 


2 


3 


+ 


T^O 


Gentile's statistics (finite p)j27] 


2 


3 


+ 


T^O 


3D Fermi paramagnet [36 


3 


3 


+ 


T^O 



Table 1: The thermodynamic curvature for several thermodynamic systems. The 
number of independent thermodynamic parameters n, spatial dimension d, sign of 
i?, and possible divergences is given. For some systems d is not determined, which 
is denoted by • • • . The sign of R is consistent with the sign convention of Weinberg 
|30j . Small \R\ means that the value of \R\ is on the order of the volume of an 
interparticle spacing or less. 
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Figure Captions 



Figure 1: An imaginary partition inside an infinite, closed universe Aq. The partition 
encloses a system A, with constant volume V. Ar is the reservoir of A, which is 
everything in Aq outside the partition. 



Figure 2: Two concentric, imaginary partitions enclosing systems Ai and A, respec- 
tively {A is part oi Ai). Ar is now everything in Aq outside of Ai. A\ samples the 
state of Aq, but A samples only the present state of Ax. 



Figure 3: The ideal gas probability P{N) of finding A^ particles in A for V such 
that (A^) = 5. P{N) is calculated with the Gaussian approximation (solid curve) 
and the exact Poison distribution (points). The Gaussian approximation is effective 
even for V close to atomic sizes. 



Figure 4: A surface on which the mean local density p equals the overall mean 
density po- Also shown is an arbitrary line intersecting the surface at the dotted 
points, and a "droplet" of linear dimension ^ equal to the mean distance between 
those intersections. 



Figure 5: For large V the density fluctuations are given by a Gaussian probability 
P(p) centered on po- If V is less than C,, fluctuations are bimodal, because A is most 
probably either inside a dense cluster or in a depleted in-between region. 



Figure 6: A two-sphere with radius a. A circle centered at the north pole is drawn 
on it. Its radius r measured along a geodesic on the two-sphere is greater than its 
radius measured in the flat three-dimensional embedding space, so its Riemannian 
curvature scalar R is negative. 
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Figure 7: The pseudosphere defined by the parametric functions x = a sin -u cos f, 
y = asinusmv, and z = a[cosM + lntan(M/2)], has constant curvature R = 2/a^. 
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